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NUMERICAL  SIMULATION  OF  A  POS8IBLE  FREEZING  AND  SHEET  FORMATION 
MECHANISM  FOR  BARIUM  CLOUD  8TRIATIONS 

1 .  Introduction 

Although  many  features  of  barium  cloud  atrlations  are  understood  (Ossakov 
1979;  Ossakow  et  al.,  1982),  there  exist  several  aspects  of  the  phenomenology 
which  as  of  yet  lack  a  satisfactory  explanation.  This  paper  addresses  two  of 
these  phenomena:  1)  the  persistence  of  strlatlons  with  scale  sizes  in  the 
kilometer  range  for  times  on  the  order  of  hours  [Prettle  et  al.t  1977,  J.A. 
Fedder  and  W.  Chestnut,  private  communication,  1980];  and  2)  the  formation  and 
persistence  of  long  thin  "sheets"  of  barium  on  the  leading  (non-structuring) 
side  of  the  cloud.  In  what  follows  we  shall  attempt  to  show  that  both  of 
these  phenomena  can  be  explained  by  a  proper  understanding  of  the  role  played 
by  the  finite  Pedersen  mobility  both  of  the  barium  ions  and  of  the  ions 
constituting  the  ambient  ionosphere.  In  essence  we  find  that  the  finite 
Pedersen  mobility  of  the  ions  vis  a  vis  the  electrons ,  which  essentially  have 
no  Pedersen  mobility  at  all,  allows  the  barium  to  "leak"  out  of  the  electron 
cloud*  Quasineutrality  in  the  electron  cloud  is  maintained  by  incoming 
ambient  ionosphere  ions  (denoted  here  by  0+,  but  these  could  also  be  N0+  or 
other  ionspheric  ion  constituents),  which  compress  up  to  the  electron  density 
as  they  replace  the  barium  in  the  cloud#  For  many  ionospheric  parameters, 
once  the  barium  has  left  the  electron  cloud,  it  is  essentially  free  of  the 
polarization  fields  associated  with  the  electron  cloud  and  hence  Immune  to 
further  structuring.  An  observer  watching  only  the  barium  might  conclude  that 
it  had  "frozen".  Further,  the  electron  cloud  itself  may  decay  if  the  ions 
coexisting  with  it  are  subject  to  fast  recombination  chemistry. 

The  above  ideas  were  arrived  at  Independently  by  Prettle  [1981]  and  by 
Fedder  (private  communication,  1981).  The  conclusions  of  the  present  paper 
differ  from  those  in  Prettle  [1981]  primarily  in  the  fact  that  we  find  that 
the  formation  of  sheets  is  directly  related  to  the  leaking  of  barium  from  the 
electron  cloud,  as  just  discussed,  and  is  totally  unrelated  to  the  formation 
of  "quadrupole"  potential  fields,  as  postulated  by  Prettle  [1981].  Rather, 
the  sheet  forms  as  the  barium  leaks  out  of  one  side  of  the  electron  cloud.  In 
the  process  it  also  leaves  the  regions  of  shielding  provided  by  the  enhanced 
Pedersen  conductivity  associated  with  the  enhanced  electron  density  (for  many 
ionospheric  parameters),  and  over  a  fairly  short  distance  enters  the 
unshielded  region  where  it  takes  on  the  ambient  c  x  B/B*  drift  velocity. 
From  the  reference  frame  of  the  center  of  the  electron  cloud,  we  are  slowly 
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injecting  barium  into  a  flow  field  transverse  to  the  injection  direction  which 
starts  at  zero  at  the  electron  cloud  center  and  rapidly  becomes  large  over  a 
very  short  distance  at  the  "edge"  of  the  electron  cloud*  The  result  is  a 
long i  low  density  string  of  barium  which  we  perceive  as  a  "sheet".  This 
concept,  as  well  as  the  meaning  of  the  term  "many  ionospheric  parameters"  used 
above,  will  be  given  in  more  detail  later  in  this  paper. 

In  section  2  of  this  paper  we  briefly  review  the  theory  and  relevant 
equations  for  the  multilevel/multispecies  ionospheric  barium  cloud  problem. 
In  section  3  we  describe  in  more  detail  the  basic  Pedersen  drift  mechanism 
proposed  by  Prettie  [1981]  and  by  Fedder  (private  communication  1981),  and 
also  present  our  proposed  mechanism  for  barium  sheet  formation.  In  section  4 
we  discuss  the  effect  of  the  Pedersen  drift  mechanism  under  various 
ionospheric  parameter  regimes  within  the  context  of  a  two  level  (two  species) 
model.  In  section  5  we  present  preliminary  numerical  simulation  results  which 
seem  to  bear  out  the  analysis,  although  one  aspect  of  the  simulations,  the 
steep  density  gradients  on  one  side  of  the  sheet,  is  not  well  understood  by  us 
at  this  time.  Finally  in  section  6  we  present  our  conclusions  and  plans  for 
future  work. 


2.  Theory 

Consider  a  fluid  plasma  species  a  imbedded  in  a  neutral  gas  and  a 
magnetic  field  B,  which  we  assume  is  aligned  along  the  cartesian  axis  z. 
Under  the  time  scales  and  ionospheric  parameters  of  interest  for  the 
ionospheric  barium  cloud  strlations  problem,  this  plasma  will  respond  to  an 
external  force  perpendicular  to  the  magnetic  field,  JF^,  *n  two  ways:  1)  by 
drifting  in  the  direction  of  the  external  force  (Pedersen  mobility);  and  2)  by 
drifting  in  a  direction  perpendicular  to  both  the  external  force  and  to  the 
magnetic  field  (Hall  mobility).  Specifically,  if  the  momentum  equation  for 
species  a  is 


+  v  *7lv  »  —  (E  +  )  “  v  (v  -  U  )  +  —  G 
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and  we  define  Che  external  force  to  be 
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Here  we  have  neglected  the  left  hand  side  of  (2.1)  since  we  are  Interested 

only  In  average  drift  velocities  over  tine  scales  that  are  long  coapared  to 

either  the  mean  time  between  collisions  or  the  gyroperiod.  Also  we  have 

ignored  all  collisions  except  those  between  species  a  and  the  neutral  gas. 

The  notation  in  the  above  equations  is  as  follows:  and  m^  are  the 

fluid  velocity,  charge  per  particle,  and  mass  per  particle  respectively  of 

species  o;  E_ is  the  electric  field,  t  is  time,  c  is  the  speed  of  light,  v  is 

the  collision  frequency  of  species  o  with  the  neutral  gas,  U„  is  the  neutral 

gas  velocity,  fla  is  the  cyclotron  frequency  of  species  a,  0^  =  Iq^B/ (m^c) I , 

and  G  is  a  term  meant  to  represent  all  of  the  other  external  forces  on 
a 

species  o  (gravity,  pressure  gradients,  etc).  The  subscript  l  refers  to  the 
components  of  the  respective  vector  which  are  perpendicular  to  _B. 

Starting  with  the  above,  one  can  derive  the  set  of  partial  differential 

equations  describing  the  barium  cloud  striatlons  problem  by  demanding  tl  it 

each  speciea  a  obey  a  continuity  equation  and  that  -  0  where  is  the 

electric  current,  which  is  to  say  that  tie  demand  electrical  quasi-neutrality 

globally  for  all  time.  The  quantity  is  defined  by  multiplying  (2.3) 

by  q  n  ,  where  n  is  the  number  density  of  species  a,  and  sunning  over  a.  If 
a  a  a 

one  further  assumes  that  ■  -  V$  where  4  is  the  scalar  electrostatic 

potential,  then  V»J[  •  0  yields  an  elliptic  equation  for  p.  This  derivation  is 
given  in  several  papers  f Lloyd  and  Haerendel,  1973;  Goldman  et  al»,  1974; 
Perkins  et  al.,  1973  for  example]  and  will  not  be  reproduced  here.  Rather,  we 
simply  describe  the  other  simplifying  assumptions  and  give  the  final 
equations. 
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We  find  that  for  ionospheric  barium  clouds  it  is  a  very  good 
approximation  to  assume  that  the  electrostatic  potential  t  is  constant  along 
magnetic  field  lines,  i.e.  $  *  $(x,y)  where  the  magnetic  field  is  aligned 
along  t  in  an  x-y-z  cartesian  coordinate  system.  Further  we  find  that  since 
the  currents  parallel  to  are  carried  primarily  by  electrons  and  since  the 
motion  of  ions  parallel  to  B^  consists  primarily  of  a  slow  diffusion  plus  a 
bulk  “falling"  of  the  cloud,  since  *  0  where  £  is  the  gravitational 
acceleration,  it  is  sufficient  to  represent  the  ions  as  an  array  of  two- 
dimensional  planes  of  plasma  perpendicular  to  B_,  each  moving  with  the  bulk 
"falling"  velocity  along  ^B,  and  hence  to  treat  numerically  only  the  transport 
of  ions  perpendicular  to  B_  within  each  layer.  For  cases  where  multiple  Ion 
species  coexist  in  one  layer,  this  layer  is  simply  treated  as  multiple  layers 
with  the  same  location  In  space.  This  treatment  is  consistent  with  our 
neglect  of  all  collisions  except  ions  with  neutrals  and  electrons  with 
neutrals.  Thus  we  use  the  phrase  "multilevel/multispecies"  to  describe  the 
model. 

We  draw  the  readers  attention  to  the  following  remarks,  whose  Importance 
shall  become  clear  as  we  progress.  The  right  hand  side  of  (2.3)  gives 
explicitly  the  Pedersen  and  Hall  components,  respectively,  of  the  drift 
velocity  of  species  a.  Looking  at  (2.4)  we  find  that  k^  *  0  for  electrons, 
since  ven/^e  m  0  (where  the  subscript  e  denotes  electrons).  Thus  ions,  for 
which  van/^a  ranges  from  0.01  to  greater  than  1.0,  are  endowed  with  a  Pedersen 
mobility  which  is  lacking  for  the  electrons:  Ions  and  electrons  do  not  move 
together.  Equally  important,  under  the  conditions  that  Is  uniform  and  that 
the  external  force  can  be  written  as  the  gradient  of  a  scalar,  conditions 
which  are  satisfied  in  most  of  the  barium  cloud  problems  we  shall  consider, 
the  Hall  component  of  v^  is  divergence  free  while  the  Pedersen  component  is 
not.  Thus  individual  fluid  elements  of  ions  may  compress  or  expand,  while  the 
electron  fluid  satisfies  the  constraint  of  incompressibility  (in  the  x-y 
plane). 

In  this  paper  we  shall  limit  our  discussion  to  a  model  consisting  of  only 
two  layers  or  species  of  ions,  as  this  will  considerably  simplify  the 
understanding  of  the  ideas  we  shall  put  forth  without,  we  believe, 
compromising  the  ideas  themselves.  We  shall  also  confine  our  external  force 
to  an  externally  Imposed  electric  field  E^  which  is  perpendicular  to  JJ,  i.e., 


E0  -  Eox  x  +  E0y7*  We  denote  by  and  n2  the  integrated  (along  B)  number 
density  of  the  ions  (and  by  quasi-neutrality,  the  electrons)  in  levels  1  and  2 
respectively*  The  equations  describing  this  system  are  then 

Bn^/at  +■  “  0  0  ■  !»2  (2*6) 

V«v  +  Wl  * B  *  W£Pi  +  £p2l  <2-7> 

B  5  -  k  «£hi  +  S*Kfc  -  v»  +  h  ttEhi  ♦  £h2KH  -  £oJ)  <2-8> 
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a 

v 
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a 

The  terms  comprising  H  in  (2.7)  are  referred  to  as  the  "Hall  terms",  and  for 
the  purposes  of  the  ideas  to  be  presented  in  the  present  paper  can  be 
considered  to  be  small.  The  quantity  in  (2.6)  is  defined  by  (2.3)-(2.5) 
where 

loti  -  q0<Io  -  w>  (2.11) 

3.  The  Pedersen  Leakage  Mechanism:  The  Simplest  Case 

As  mentioned  in  the  previous  section  we  shall  consider  only  a  model 
consisting  of  two  ion  species.  Since  we  are  neglecting  all  collisions  except 
those  between  ions  and  neutrals  (recall  we  have  set  ven/^e  *  0) >  «nd  since  we 
are  assuming  that  the  only  external  forces  acting  are  electric  fields 
perpendicular  to  B,  which  are  constant  along  any  given  magnetic  field  line, 
one  can  conceptually  think  of  these  layers  of  plasma  as  coexisting  at  the  same 
plane  in  space,  as  long  as  we  enforce  the  caveat  that  the  two  ion 
species  o  may  "see**  different  neutral  densities,  which  will  be  reflected  in 
the  value  of  vaQ.  Hence  all  of  our  diagrams  of  the  physical  processes  can  be 


of  a  single  two-dimensional  plane* 

Furthermore,  we  shall  specifically  take  the  two  ion  species  to  be  Ba*  and 
the  ambient  ionosphere  ion  0*,  with  the  understanding  that  by  "0*M  we  might 
also  mean  NO*  or  any  of  the  other  ambient  ionospheric  ions.  Thus  we  shall 
henceforth  use  the  subscripts  b  and  o  to  denote  Ba*  and  0*,  rather  than  the 
subscripts  1  and  2  used  previously.  Within  the  context  of  this  model,  it  is 
clear  that 


nb  + 


(3.1) 


where  ne  denotes  the  Integrated  electron  density.  Hence  one  can  substitute  a 
continuity  equation  for  nfi 


3ne/3t  +  V(nevel>  *  0  (3*2) 

for  one  of  the  two  equations  In  (2.6),  since  the  substituted  density  can 
always  be  recovered  from  (3.1).  For  the  remainder  of  this  paper  we  shall  use 
this  fact  and  display  profiles  of  barium  density  and  electron  density*  For  the 
purposes  of  this  section  only,  ve  impose  one  last  assumption: 


»„/a. 

on  o 


(3*3) 


Thus  the  two  ion  species  are  identical  in  terms  of  their  electrical  properties 
vis  a  vis  (2.7)  and  in  terms  of  their  response  to  an  external  force  (2.3).  In 
fact,  if  we  were  interested  only  in  electron  densities  and  not  in  the  actual 
Identities  of  the  Ion  species  making  up  the  "electron  cloud**,  it  would  suffice 
to  sum  the  two  ion  layers  and  treat  them  as  simply  one  layer*  This  is 
precisely  the  implicit  assumption  made  in  the  "one  level  model"  [McDonald  et 
al*,  1980]*  Instead  we  shall  choose  to  follow  each  ion  species  separately,  as 
this  will  have  interesting  consequences. 

In  Fig.  1  we  show  a  simplified  picture  of  an  ionized  barium  cloud  which 
has  been  placed  in  a  uniform  ionosphere  consisting  of  0*  ions*  Isodensity 
contours  of  electrons  are  shown  as  solid  lines,  while  those  for  barium  are 
shown  as  dashed  lines.  At  t  •  0  the  two  sets  of  contours  are  identical.  We 
Impose  an  ambient  electric  field  directed  to  the  right*  The 
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electron/barium  cloud  will  partially  shield  Itself  from  this  field,  and 

experience  an  average  reduced  field  E' ,  still  directed  to  the  right  In  Fig. 

1.  Referring  to  (2.3)-(2.5)  and  recalling  that  v  /fl  «  0  we  find  that  the 

en  e 

electron  cloud  simply  drifts  downward  at  velocity  cE'/B •  If  we  assume  that 
-2 

Vbn^b  i8  ~  »  then  we  see  that  the  barium  cloud  too  drifts  downward  at  a 

velocity  very  close  to  cE'/B «  However  the  barium  also  experiences  a  Pedersen 

drift  to  the  right  of  value  approximately  cE'/B .  If  the  clouds 

Initially  have  diameter  D  and  remain  undeformed  the  barium  cloud  will 

completely  exit  the  electron  cloud  in  a  time  t  -  (fl./v.  )DB/(cE"),  as  depicted 

D  bn 

on  the  left  side  of  Fig.  1.  Of  course  by  quasi-neutrality  the  ions  in  the 
electron  cloud  would  then  consist  of  0+  ions  which  had  been  compressed  up  to 
the  required  density  (the  electron  density)  as  they  entered  the  electron  cloud 
from  the  left*  By  the  same  principle,  the  barium  must  expand  as  it  leaves  the 
electron  cloud  to  the  right,  showing  that  this  motion  cannot  truly  take  place 
without  deformation*  Also  one  might  expect  to  see  the  normal  electron  cloud 
steepening  and  structuring  during  this  time  interval.  In  addition  as  the 
barium  exits  the  electron  cloud  to  the  right  it  leaves  a  region  where  the 
electric  field  has  been  reduced  by  shielding  to  E'  and  enters  a  region  where  E 
takes  on  the  larger  value  EQ.  Accordingly  as  it  leaves  it  gets  caught  up  in 
this  higher  velocity  field  and  forms  a  sheet  of  lower  density  barium  which 
leads  the  cloud.  This  somewhat  more  realistic  picture  is  depicted  on  the 
right  side  of  Fig.  1. 

The  above  picture  describes  In  qualitative  terms  the  essence  of  this 
paper.  We  will  now  expand  upon  this  idea  in  more  quantitative  terms  by 
simplifying  the  picture  even  further  so  that  quantitative  results  can  be 
obtained  and  that  a  more  precise  picture  may  be  obtained  of  the  compression 
and  expansion  processes  proceeding  in  the  ion  fluids. 

In  Fig.  2a  we  again  show  a  profile,  n^  vs  x,  of  an  ionized  barium  cloud 
which  has  been  placed  in  a  uniform  ionosphere  consisting  of  0+  ions,  again 
subject  to  an  externally  imposed  electric  field  E^  pointing  to  the  right. 
This  time,  however,  we  choose  the  cloud  and  ionosphere  to  be  Infinite  In  the  y 
direction,  the  so-called  slab  geometry.  We  denote  the  cloud  thickness  by  L 
and  assume  for  the  sake  of  example  that  its  Integrated  density  is  exactly 
twice  that  of  the  background  0+  ionosphere.  The  barium  fluid  is  shown  cross- 
hatched,  while  the  0+  fluid  is  shown  without  cross-hatching.  Again  we  are 
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BARIUM  CLOUD 
ELECTRON  CLOUD 


SIMPLIFIED 


MORE  REALISTIC 


(Q/yjP 

(cE'/B) 


ELECTRON  BARIUM 
CLOUD  CLOUD 


ELECTRON 
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Fig.  1  —  Schematic  depiction  of  the  evolution  of  an  ionized  barium  cloud  subject 
to  an  externally  imposed  rightward  directed  electric  field  Eq.  The  cloud  will 
partially  shield  itself,  yielding  a  net  field  E'  inside  the  cloud.  Rightward  Pedersen 
drifts  then  allow  the  barium  to  separate  from  the  electron  cloud  in  a  time  t<,. 
A  simplified  picture  is  shown  on  the  left,  while  a  more  realistic  evolution  is 
depicted  on  the  right. 
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Then  the 
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_2 

assuming  that  v.  /fl,  ■  v  /ft  and  that  both  quantities  are  ~  10 
bn  d  on  o 

potential  equation  (2*7)  has  the  simple  solution 


E  -  E 

o 

-  V3 


outside  the  electron  cloud 
inside  the  electron  cloud 


(3.4) 


If  is  taken  to  be  upward  in  Fig.  2  then  both  ion  species  and  the  electrons 
are  all  drifting  toward  the  reader  with  velocity  very  close  to  cE/B.  However 
both  the  ions  have  in  addition  a  Pedersen  drift  to  the  right  given  by 


»p  ■  (”bn/VcE/B  '3-5> 

where  vp  is  the  rightward  Pedersen  velocity  in  this  example.  Noting  the 
discontinuous  changes  in  E  that  occur  at  the  edges  of  the  electron  cloud,  we 
see  that  ion  compression  is  taking  place  at  the  left  edge  of  the  electron 
cloud  while  expansion  is  taking  place  at  the  right  edge.  In  fact  in  this 
example ,  the  ions  experience  a  compression  of  exactly  a  factor  of  3  as  they 
cross  the  left  edge  of  the  electron  cloud,  and  an  expansion  of  exactly  a 
factor  of  3  as  they  exit  the  right  edge.  Of  course  this  is  exactly  what  is 
required  by  the  condition  of  quasi-neutrality  (7*^  ■  0).  One  can  also 
calculate  the  time  tQ  it  takes  for  the  barium  to  completely  exit  the  electron 
cloud  in  a  manner  analogous  to  that  used  previously  for  Fig.  1.  One  obtains 


t 

o 


v.  cE 
f  bn\  _ o 

J  38 


(3.6) 


where  L  Is  the  width  of  the  cloud.  In  Fig.  2b  we  show  the  tine  evolution  of 

the  profile  depicted  In  Fig.  2a  for  various  units  of  tine  In  terns  of  tQ. 

Note  that  after  tine  tQ:  1)  the  "cloud"  consists  totally  of  0*  Ions;  2)  the 

barium  extends  over  a  region  three  tines  Its  original  extent  but  down  by  a 

factor  of  three  in  density;  3)  the  barlun  Is  conpletely  renoved  from 

gradients  In  na  and  hence  £  and  therefore  would  be  renoved  from  any  further 
*  P 

structuring  taking  place  In  the  analogous  situation  for  a  two-dinenslonal 
cloud,  as  In  Fig.  1.  It  is  postulated  that  this  may  very  well  be  what  is 

known  as  "freezing".  It  should  also  be  noted  that  if  the  0*  cloud  is  subject 


COMPRESSION  EXPANSION 


t0  =  TIME  FOR  Ba+  IONS  TO  LEAVE  SLAB 
L 


Fig.  2(a)  —  Schematic  of  a  barium  cloud  in  an  ionosphere  consisting  of  0+ 
ions,  similar  to  Fig.  1  except  that  now  we  are  in  a  “slab”  geometry  (plasma  has 
infinite  extent  in  y  direction).  The  barium  ions  are  shown  cross-hatched  while  the 
oxygen  ions  are  shown  without  cross-hatching.  The  shielding  by  the  cloud  causes 
a  compression  of  ions  to  take  place  at  the  left  edge  and  an  expansion  to  take  place 
at  the  right  edge.  The  envelope  n,,  +  nb  depicts  the  electron  density  i^,  which 
cannot  change  in  time. 
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to  short  chemical  recombination  times,  this  mechanism  may  very  well  result  in 
the  decay  of  the  electron  density  cloud  as  well  as  the  freezing  of  the  barium 
cloud*  This  recombination  effect  is  discussed  in  more  detail  in  Prettle 
[1981]. 


We  are  now  in  a  position  to  understand  how  it  is  that  the  barium  that 

leaks  from  the  electron  cloud  via  the  above  mechanism  appears  in  the  form  of 

sheets  of  barium  leading  the  electron  density  cloud.  Referring  to  Fig.  2a,  we 

note  that  as  the  barium  exits  the  right  edge  of  the  electron  cloud,  tt 

experiences  a  factor  of  three  increase  in  the  rightward-directed  electric 

field.  It  is  the  change  in  Pedersen  drift  associated  with  this  which  is 

responsible  for  the  expansion  of  the  ions  as  they  exit  the  electron  cloud. 

However  it  is  precisely  this  same  change  in  E  which  will  produce  a  three-fold 

increase  in  the  cE/B  velocity  of  the  ions  toward  the  reader.  Now  consider  the 

slice  of  barium  located  at  y  -  0  at  t  -  0.  The  right  hand  edge  of  this  slice 

will  immediately  start  moving  toward  the  viewer  at  a  velocity  cEQ/B  and  to  the 

right  at  velocity  (v  /ft  )  cE  / B.  The  left  hand  edge  of  this  slice  will  move 

on  d  o 

at  precisely  these  same  velocities  but  reduced  by  a  factor  of  three.  At  time 

tQ  when  the  barium  has  completely  exited  the  electron  cloud  of  width  L,  this 

slice  will  have  Its  right  hand  and  left  hand  coordinates  located  at 

(x  -  3L,  y  -  -  3  ft.v  *  L)  and  (x  -  0,  y  ■  -  ft,v  *  L)  respectively,  where  we 

b  bn  b  bn 

have  defined  x  -  0  to  be  the  right  edge  of  the  electron  cloud.  The  barium 
slice  is  now  31  wide  in  the  x  direction  and  2L  long  In  the  y 

direction.  For  ft^v^  -  20,  we  get  a  barium  slice  which  is  13  times  larger  in 
y  than  it  is  in  x.  Referring  to  Fig.  I  we  see  that  our  postulated  slice  in 

Fig.  2  is  analogous  to  the  barium  cloud  itself  in  Fig.  1,  so  that  again  we 

expect  the  barium  which  has  leaked  out  of  the  electron  cloud  to  be  much  longer 
in  the  y  direction  than  in  the  x  direction,  a  geometry  we  could  easily  term  a 
"sheet” .  The  process  is  quite  similar  to  what  would  happen  if  one  were  to 
eject  a  column  of  smoke  sidewards  from  a  moving  automobile:  the  trail  of 
smoke  would  be  much  larger  along  the  direction  of  travel  than  tranverse  to 
it.  From  the  frame  of  reference  of  the  electron  cloud  in  Fig.  1,  we  are 
injecting  a  column  of  barium  "smoke”  transverse  to  our  direction  of  travel 
relative  to  the  ambient  plasma.  Hence  the  "trails"  of  barium. 


Prettle  [1981]  has  suggested  that  the  barium  sheets  are  the  result  of 
image  formation  and  the  resultant  "quadrupole"  electric  fields.  While  we 


H 


would  be  the  last  to  suggest  that  Image  formation  is  not  important  in  barium 
cloud  striations  [Scannapieco  et  al.,  1976;  Ossakow  et  al.t  1980],  we  believe 
the  above  shows  the  explanation  to  be  quite  a  bit  simpler.  Numerical 
simulations  to  be  presented  in  section  5  also  bear  this  out. 


4.  The  Pedersen  Leakage  Mechanism:  More  General  Ionospheric  Parameters 

In  the  last  section  we  confined  ourselves  to  the  case  of  two  ion  species, 

each  of  which  had  the  same  value  of  v  /ft  •  In  this  section  we  relax  this 

an  a 

last  constraint  for  the  purposes  of  a  short  qualitative  discussion.  A  more 
detailed  discussion  of  this  topic,  as  well  as  a  treatment  of  the  case  of  more 
than  two  levels  will  be  given  in  a  future  paper. 

To  illustrate  the  principles  involved,  we  shall  again  use  the  one¬ 
dimensional  slab  depicted  in  Fig.  2,  with  the  understanding  that  there  is  no 
longer  a  one-to-one  correspondence  between  electron  density  ne  and  total 
Pedersen  conductivity  =  E^  +  Ep q •  Hence  we  shall  interpret  Fig.  2  to 
depict  the  ne  profile,  with  the  Ep  profile  depending  on  the  ion  species 
Involved,  and  which  is,  as  we  shall  see,  time  dependent. 


The  principle  to  have  in  mind  when  analyzing  the  general  case  is  that  the 
electrons  cannot  move  to  the  left  or  right,  since  they  have  no  Pedersen 
mobility.  Hence  the  electron  density  profile  ne  is  invariant  in  time  as 
depicted  in  Fig.  2b.  Since  n^  +  nQ  -  ne,  the  Ba+  and  0*  ions  are  constrained 
to  move  In  such  a  way  as  to  maintain  this  invariance  in  n0.  In  addition,  the 
velocity  of  each  ion  species  must  be  proportional  to  the  electric  field  and  to 
the  appropriate  van/fla-  Furthermore,  the  electric  field  must  be  inversely 
proportional  to  Ep. 

The  problem  posed  above  for  the  general  case  *  von^o  can 

solved  exactly.  We  shall  not  give  the  proof  here,  but  simply  state  the 

remarkably  simple  result:  The  time  evolution  depicted  in  Fig.  2b  for  the 

special  case  v.  /ft.  •  v  /ft  is  in  fact  the  general  solution  for 
bn  b  on  o 

Vbn^^b  *  von^o*  except  for  a  change  In  temporal  scaling.  Specifically, 


n'  «  f  n 
o  o 


(4.1) 

(4.2) 
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O//  o  ,  ox 
=  no/(no  +  nb) 

(4.7) 

M° 

i  (n°v  fl_l  +  nfv,  a”1)/(n°v  fl_1) 

o  on  o  d  bn  b  '  '  o  on  o 

(4.8) 

where  the  prime  superscript  on  a  denotes  quantities  evaluated  within  the 
expanded  barium  cloud  at  time  t  >  t0  shown  in  Fig.  2b,  the  double  prime 
superscript  denotes  quantities  within  the  electron  cloud  at  t  >  tQ  in  the  same 
figure,  the  zero  superscript  denotes  quantities  at  time  t  «  0,  and  t0  is,  as 
before,  the  time  taken  for  the  barium  cloud  to  leave  the  electron  cloud.  Note 
that  this  says  that  the  barium  will  always  leave  the  electron  cloud,  but  since 
tQ  is  proportional  to  t*ie  process  may  take  an  extremely  long  time 
if  M°  is  large  or  v^n/^  *s  small.  Just  as  important,  it  says  that  the 
Pedersen  conductivity  ratio  M  of  the  barium  and  electron  clouds  to  the  ambient 
ionosphere  will  change  in  the  general  case: 


f  M° 

(4.9) 

f1 

(4.10) 

where  the  superscripts  have  the  same  meaning  as  above.  Since  0  <  f  <  1,  and 
in  some  cases  may  be  close  to  zero,  the  barium  cloud  will  always  suffer  a 
decrease  in  M  value,  and  in  fact  it  is  possible  to  have  M'  <  l,  in  which  case 
the  barium  cloud  would  actually  represent  a  Pedersen  conductivity  depletion. 
As  for  the  electron  cloud,  its  M  value  may  decrease  or  Increase  or  stay  the 
same,  depending  on  f.  Sufficiently  small  f  will  result  in  large  Increases  in 
M  for  the  electron  cloud,  and  is  the  basis  for  the  freezing  mechanism  for  the 
electron  cloud  proposed  by  Ossakow  et  al.  [1981],  since  a  large  H  has  been 
shown  to  be  associated  with  a  resistance  to  further  structuring  on  the  part  of 
the  electron  cloud  [McDonald  et  al.,  1981]. 
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The  above  examples,  and  the  parameter  study  suggested  by  (4*1)  to  (4*10) 
provide  us  with  a  vivid  indication  of  the  richness  of  the  physics  in  the 
general  case  of  vki/flb  *  vOI/^o*  Thi8  ric*me8S  *8  further  compounded  when  we 
model  the  physics  using  more  than  two  layers  or  species  of  ions.  We  are 
presently  undertaking  these  parameter  studies  to  determine  these  multilayer 
effects  on  the  morphology  of  barium  cloud  striatlons.  However,  the  point  of 
this  section  is  not  to  show  how  different  and  interesting  this  regime  is,  but 
rather  to  show  that  there  are  certain  aspects  of  the  physics,  namely  the 
leakage  of  barium  from  the  electron  cloud  and  its  subsequent  reduction  in 
density  and  Pedersen  conductivity  by  a  factor  f  (Eq.  (4.1)  and  (4.10)),  that 
are  independent  of  whether  v|)n/^)  *  Von^o* 

5.  Numerical  Simulations 

In  this  section  we  return  to  the  case  considered  in  section  3,  that  of 
v,  /ft.  *  v  /ft  •  There  are  several  reasons  for  this  :  First,  we  have  shown 
in  the  previous  section  that  the  basic  Pedersen  mobility  leakage  mechanism 
proposed  here  is  active  independent  of  this  assumption.  Second,  this  is  the 
implicit  assumption  made  in  the  "one-level  approximation",  and  it  will  be 
interesting  to  see  what  effects  are  missed  when  one  only  follows  electrons. 
Finally  this  is  the  case  for  which  "image  effects"  are  absent,  and  It  becomes 
impossible  to  argue,  as  does  Prettle  [1981],  that  image-induced  quadrupole 

fields  affect  the  formation  of  barium  sheets.  As  the  numerical  simulations 

will  show,  well-defined  sheets  of  barium  are  produced  without  any  evidence  of 
quadrupole  fields. 

The  simulation  was  performed  on  a  stretched  grid  with  43  and  120  grid 
points  in  the  x  and  y  directions  respectively.  The  central  31  by  100  portion 
of  the  mesh  was  uniformly  gridded  with  Ax  *  Ay  ■  0.258  km.  The  four 

boundaries  of  the  grid  were  placed  well  away  from  the  cloud  (utilizing  the 
stretched  mesh)  to  ensure  no  influence  of  the  boundaries  on  the  evolution  of 
the  cloud.  Boundary  conditions  on  all  plasma  species  were  that  the  normal 
derivative  vanished.  Boundary  conditions  on  $  were  Dirichlet  ($  »  0)  on  the 
left  and  right  boundaries,  and  Neumann  (d^/dy  •  0)  at  the  top  and  bottom 
boundaries.  Two  ion  species  were  used,  Ba+  and  0+,  with  a  value 

of  van/fla  “  0.06  assigned  to  each.  The  quantity  nQ  was  initially  assigned  a 
normalized  value  of  1.0,  and  the  Initial  barium  cloud  is  given  by 
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(5.1) 


nb  -  4.0  exp(-  r4/L4)(l+P(x)) 

r2  -  (x-xQ)2  +  (y-y0)2  (5.2) 

L  -  2.5  kn  (5.3) 

P(x)  -  0.05  cos(20(x-xo)/L)  (5.4) 

where  the  point  (x0,yQ)  specifies  the  Initial  center  of  the  cloud.  The 

magnetic  field  Is  taken  to  be  0.5  gauss  directed  along  the  positive  * 
axis.  The  ambient  electric  field  ^  Is  directed  along  the  positive  *  axis 
with  a  magnitude  such  that  c^  x  _B/B2  is  100  m/sec  directed  along  the  negative 
9  axis.  The  numerical  resolution  of  this  cloud  is  lntenttonally  low,  in  order 
to  allow  us  to  run  the  calculation  to  late  times  at  moderate  cost.  The 

Inherent  numerical  dissipation  of  this  low  resolution  has  the  effect  of 
Inhibiting  bifurcation  of  the  cloud  but,  as  we  shall  see,  this  Is  actually  an 
advantage  since  It  allows  us  to  see  the  leakage  of  the  barium  from  the 
electron  cloud  In  a  "clean”  calculation,  without  the  complicating  effects  of 
cloud  bifurcation. 

Figure  3  shows  isodensity  contours  of  ne  (cross-hatched,  denoted  "Density 
1")  and  nb  (denoted  "Density  2")  at  various  times  during  the  calculation.  The 

slight  asymmetry  in  the  electron  density  plots  is  due  to  the  Hall  terms  In  the 

potential  equation.  All  calculations  are  performed  in  a  reference  frame 
moving  with  the  electron  cloud  center  of  mass.  Note  that  at  time  t  -  0  sec 
the  barium  and  electron  contours  are  centered  on  one  another,  but  that  as  time 
progresses  the  barium  shows  its  inevitable  drift  to  the  right,  "leaking"  out 
of  the  electron  cloud.  At  the  latest  time  shown,  1096  sec,  virtually  all  of 
the  barium  has  left  the  electron  cloud  and  formed  a  long,  low  density  sheet  of 
barium  plasma  to  the  right  of  and  leading  the  electron  cloud.  One  feature  of 
the  sheet  that  is  quite  striking  is  the  much  larger  transverse  gradients  of  nb 
on  the  left  side  of  the  sheet  than  on  the  right  side.  There  exists  the 
possibility  of  optical  or  radar  measurements  of  this  barium  cloud  interpreting 
this  steep  gradient  as  "frozen  structure". 
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Fig.  3  (Cont’d)  —  Evolution  of  the  barium  cloud  described  in  the  text  computed  via  numerical 
simulation.  Electron  density  contours  (cross-hatched)  are  shown  on  the  left,  while  contours  of 
barium  density  are  shown  on  the  right,  (c)  658  seconds  and  (d)  1096  seconds. 


A  possible  explanation  of  this  feature  is  as  follows.  The  barium 
comprising  the  right  edge  of  the  sheet  originated  near  the  right  edge  of  the 
original  barium  cloud  and  thus  was  subject  to  some  expansion  as  it  left  the 
electron  cloud.  Hence  the  right  side  of  the  sheet  can  be  expected  to  have 
gradients  as  shallow  or  shallower  than  that  of  the  original  cloud.  The  barium 
comprising  the  left  edge  of  the  sheet  originated  on  the  left  side  of  the 
original  barium  cloud  and  was  therefore  compressed  in  the  x  direction  at 
first,  resulting  in  a  steepening  of  contours  in  the  x  direction.  Subsequently 
the  barium  was  expanded  but  this  time  in  the  y  direction,  a  process  which  has 
no  effect  at  all  on  the  relative  steepness  of  gradients  in  the  x  direction. 
Thus  the  steepness  of  the  contours  in  the  x  direction  is  preserved.  Fig.  4 
we  show  contours  of  constant  electrostatic  potential  4  at  the  same  times  as 
those  used  in  Fig.  3.  Notice  the  absence  of  any  structure  which  might  look 
even  remotely  like  a  "quadrupole". 

It  should  be  pointed  out  that  while  the  density  of  the  barium  comprising 
the  sheets  may  be  somewhat  low,  it  will  appear  bright  optically  to  any 
observer  whose  line  of  sight  is  approximately  parallel  to  the  sheet 
orientation  since  the  barium  is  optically  thin  and  the  integration  path 
long.  For  a  barium  cloud  driven  by  neutral  winds,  this  will  be  the  case  for 
an  observer  anywhere  near  the  initial  cloud  release  site  watching  the  cloud 
drift  away  at  late  times. 

6.  Conclusions  and  Future  Work 

We  believe  that  we  have  shown  the  Pedersen  drift  mechanism  proposed  by 
Prettle  [1981]  and  by  Fedder  (private  communication,  1981)  to  be  a  viable 
candidate  for  barium  cloud  freezing.  In  addition  we  have  shown  that  this 
mechanism  can  produce  sheets  of  barium  at  late  times  which  lead  and  are 
displaced  from  the  electron  cloud.  We  do  not  suggest,  however,  .hat  this  is 
the  only  possible  mechanism.  Indeed,  we  have  proposed  alternative  mechanisms 
ourselves,  and  are  presently  considering  still  others.  Nonetheless  the 
mechanism  looks  sufficiently  promising  to  merit  further  study.  Certainly 
higher  resolution  numerical  studies  are  in  order,  along  with  parameter  studies 
of  the  effects  of  the  relative  Pedersen  mobilities  of  the  barium  and 
ionospheric  constituents.  We  hope  to  report  on  these  in  the  near  future. 


*  calculation  shown  in  Fig.  3.  Positive  values  of  ^  are 
contoured  as  dashed  lines.  Note  the  absence  of  any 
seconds  and  (b)  contours  at  658  and  1096  seconds. 
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SAN  DIEGO,  CA  92138 
01  CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20707 

01 CY  ATTN  DOCUMENT  LIBRARIAN 
01 CY  ATTN  THOMAS  POTEMRA 
01 CY  ATTN  JOHN  DASSOULAS 
01 CY  ATTN  DR.  DONALD  J.  WILLIAMS 


KAMAN  SCIENCES  CORF 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAHER  QQ) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DAS  LAC 
01 CY  ATTN  TIM  STEPHANS 
01  CY  ATTN  WARREN  S.  KNAPP 
01  CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
01CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01CY  ATTN  DEPT  60-12 
01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 
01CY  ATTN  MARTIN  WALT  DEPT  S2-12 
01CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01CY  ATTN  RICHARD  G.  JO*  SON  DEPT  52-12 
01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  &  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  &  ENGR.  CTR. 

4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALE  H.  DIVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 
OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 


MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01 CY  ATTN  J.  HOULE 
01 CY  ATTN  GEORGE  MR0Z 
01CY  ATTN  W.  OLSON 
01  CY  ATTN  R.H.  HALPRIN 
01 CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01 CY  ATTN  W.F.  CREVIER 
01 CY  ATTN  STEVEN  L.  GUTSCHE 
01  CY  ATTN  D.  SAP  PEN  FI  ELD 
01  CY  ATTN  R.  BOGUSCH 
01 CY  ATTN  R.  HENDRICK 
01 CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01 CY  ATTN  M.  SCHEIBE 
01CY  ATTN  CONRAD  L.  LCNGMIRE 
01 CY  ATTN  WARREN  A.  SCHLUETER 

NITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01  CY  ATTN  JOHN  MORGAN  STERN 
01  CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01 CY  ATTN  W.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 
01CY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01 CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01 CY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  HA  98009 
01CY  ATTN  E.J.  PREMOUH 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

RAD  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WIGHT,  JR. 
01CY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  KA  01776 
OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LIN  SON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUCALL 

SCIENCE  APPLICATIONS,  INC 
1710  COODRIDCE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 
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SRI  L'JTERNATIONAL 
333  RAVEN SWOOD  AVENUE 
MENLO  PARK,  CA  94025 


OICY 

ATTN 

DONALD  NEIL  SON 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

L.L.  COBB 

OICY 

ATTN 

DAVID  A*  JOHNSON 

OICY 

ATTN 

WALTER  G.  CHESNUT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

M.  BARON 

OICY 

ATTN 

RAY  L.  LEADABRAND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

J.  PETERSON 

OICY 

ATTN 

R.  HAKE,  JR. 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
I  DE  ANGELO  DRIVE 
BEDFORD,  MA  01730 
OICY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ODE  COPY  TO  EACH  OP  THE 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
DR.  P.  MANGE  -  CODE  4101 
DR.  R.  MEIER  -  CODE  4141 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  GOODMAN  -  CODE  4180 
DR.  R.  RODRIGUEZ  -  CODE  4187 
CODE  2628  -  20CY 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


FOLLOWING  PEOPLE: 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAEDA 
DR.  R.F.  BENSON 

NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUCH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 
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PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR*  F.  PERKINS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  FRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  1 5213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 


U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 


UNIVERSITY  OF  CALIFORNIA 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

J-10,  MS-664 


LOS  ALAMOS,  NM 
M.  PONGRATZ 
D.  SIMONS 
G.  BARASCH 
L.  DUNCAN 
P.  BERNHARDT 
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UNIVERSITY  OF  CALIFORNIA, 
LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F. V.  CORONITI 
DR.  C.  KENNEL 
DR.  A. Y.  WONG 


UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 


